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Abstract 


The main objective of this article is to introduce and to show the 
implementation of a novel two-stage procedure to efficiently estimate the 
level of scale resolution possible for a given flow on a given grid for 
Partial Averaged Navier-Stokes (PANS) and other hybrid models. It has 
been found that the prescribed scale resolution can play a major role in 
obtaining accurate flow solutions. The first step is to solve the unsteady 
or steady Reynolds Averaged Navier-Stokes (URANS/RANS) equations. 
From this preprocessing step , the turbulence length-scale field is 
obtained. This is then used to compute the characteristic length-scale 
ratio between the turbulence scale and the grid spacing. Based on this 
ratio, we can assess the finest scale resolution that a given grid for a 
given flow can support. Along with other additional criteria, we are 
able to analytically identify the appropriate hybrid solver resolution for 
different regions of the flow. This procedure removes the grid 
dependency issue that affects the results produced by different hybrid 
procedures in solving unsteady flows. The formulation, implementation 
methodology, and validation example are presented. We implemented 
this capability in a production Computational Fluid Dynamics (CFD) 
code, PAB3D, for the simulation of unsteady flows. 

Introduction 


The limited capability of the Reynolds Averaged Navier-Stokes (RANS) approach combined with 
eddy- viscosity turbulence models to simulate unsteady and complex flows has been well known for some 
time. The RANS’ assumption is that most of the energy is modeled through the turbulence transport 
equations and is resolved in the grid. RANS also overpredicts the eddy viscosity, which results in 
excessive damping of unsteady motion. Consequently, the eddy viscosity attains unphysically large 
values due to unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved 
flow field. One of the approaches used to overcome this problem is to provide a mechanism for the RANS 
equations to resolve the largest scales of motion. Among several methods, the Detached Eddy 
Simulations (DES) [1], the hybrid Large Eddy Simulation (LES) [2-3], the Limited Numerical Scheme 
(LNS) [4] and the Partial Averaged Navier-Stokes (PANS) [5] provide the mechanism best needed to 
satisfy this requirement. 

PAB3D is a structured, multiblock, parallel, implicit, finite-volume solver of the three-dimensional 
RANS equations. Investigations in the area of unsteady flow control for propulsion applications have led 
to an increased interest in upgrading PAB3D’s [6-8] time-accurate capabilities. Advanced turbulence 
models are available in the code, which are widely used in internal and external flow applications by 
NASA and by the US aerospace industry. 

In an attempt to increase the fidelity and accuracy of the PAB3D code, a hybrid turbulence model 
RANS/LES [2-3] has been added. Another new feature to PAB3D is the addition of the PANS method, 
which was suggested by Girimaji [4]. The PANS model was developed to overcome the grid dependency 
associated with the use of the hybrid RANS/LES models. The primary objective of this paper is to 
develop a procedure which can help in the efficient implementation of PANS and other hybrid methods. 



We implemented this capability in a production Computational Fluid Dynamics (CFD) code, PAB3D, for 
the simulation of unsteady flows. The formulation, implementation methodology, and validation example 
are presented in support of the enhanced PAB3D's time-accurate and turbulence modeling capabilities. 
The new features were utilized to compute the flow field past a stationary cylinder at Re=50,000. We 
provide three-dimensional RANS and PANS simulations for this flow. Predictions of the Strouhal number 
St, the drag coefficient C D , and the lift coefficient C L are presented. We will only consider the case in 
which the flow at the separation point is laminar. This preliminary result opens the way to the potential 
capability of PANS in simulating the unsteady flow phenomena for complex aerodynamic configurations. 

The organization of the paper will be as follows: 

1) The governing equations of the RANS and PANS formulations will be presented and discussed in 
detail, and 

2) Computational results from RANS and PANS for a flow past a stationary cylinder will be 
presented, examined, and compared to experimental data. Flow around a cylinder is considered as 
a standard test case for the hybrid turbulence model [9-11], because it poses a basic flow physics 
problem and is inherently unsteady. 

Approaches 


The governing equations of the time-averaged formulation include the conservation equations for 
mass, momentum, and energy, and the equation of state. In the present study, the perfect gas law is 
chosen to represent the air properties, and the eddy-viscosity concept is used to model the Reynolds 
stresses. The mass, momentum, and energy conservation equations of the time-averaged equations can be 
written in a conservative form as follows: 
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In the case of Reynolds Averaged Navier-Stokes (RANS) equations, a Standard Turbulence Model 
(STM) such as the two-equation (ks) turbulence model is used: 
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The turbulent stress components are formulated as: 
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For the purpose of this paper, we will define RANS turbulent viscosity as 

(4) 
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The PANS model [5] was developed to overcome the grid dependency associated with the use of other 
Hybrid Turbulence Models (HTM). In its original form, PANS [5] replaces the two-equation turbulence 
model by solving for the unresolved kinetic energy k u and the dissipation e u . The k u equation is identical 
to the original k equation. In the 8 equation, the following coefficients are used to change the two- 
equation model to the Hybrid Turbulence Model (HTM), which becomes known as the PANS 
formulation through the following changes: 

A-2 = fkipe 2 ” Ql) + Ql ^ 

°k = fk°k and V = fk°e 
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The original formulation uses a constant value for the unresolved kinetic energy parameter (f k ). The users 
will select a value for this parameter and refine the grid until the flow solution converges toward a 
solution target. This could be very time consuming for resolving complex three-dimensional flows. 

Single-Stage Approach (PANS-vl) 

Elmiligui et al. [12] introduced an attempt to use a variable f k with time and space instead of using the 
constant value. In reference 12, the following equation is used to compute f k as: 

/;={!.+ tanh[2;r( A - 0.5)]} /2. (6) 


In this case, turbulent length scale is defined as: 



and the characteristic length scale ration as : 

A = f , A = maxCA^A^A.) 

A x ,A and A^are the grid cell distances in the x, y, and z directions, respectively. The function in 

equation (6) has the characteristic to be equal to 1.0 in the viscous sub-layer, as the unresolved 
characteristic ratio tends to be of very small value. Also, the value of this function is restricted to 1, in 
case the RANS turbulent viscosity becomes smaller than the LES viscosity. We define LES viscosity as: 


v EES = 0.084A VE® and k LES = f k k (7) 


Two-Stage Approach 

Some of the drawbacks of such a function used in PANS-vl are that it varies with time and space 
and that it could be affected by grid resolution. In this section, we use another approach to construct a 
function which varies only with space. 

One of the major deficiencies associated with the heretofore published use of hybrid schemes such as 
DES and PANS is that there is no clear identification of the different flow regions as shown in Figure 1. 
These regions need to be clearly defined as RANS regions and DES/PANS regions in order to achieve 
complete simulation, independent of grid resolution. Several researchers [2, 3, 12, 13 and 14] observed 
that in most cases, using hybrid methods, the use of fine grid might result in incorrect simulations. We 
believe that this is caused by the fact that the DES/PANS regions encroach into the RANS regions and 
destroy the ability of the hybrid methods to produce correct simulation. It should be understood that 
hybrid methods rely on the RANS regions in order to correctly develop the other flow regions. 

It is quite difficult to develop a method that isolates the RANS regions at all times within a hybrid 
simulation. One approach is to pre-design the grid as done by Travin et al. [13] in the case of the flow 
around a stationary cylinder. This can be accomplished for a very simple flow, where we can easily 
identify the different flow regions. In the present paper, we introduce a more general approach to define 


4 



these flow regions. This approach uses a two-stage procedure to define f k . In the first stage, we solve the 
problem using Standard Turbulence Models (STM) such as a linear two-equation model. We use the 
results of this simulation to identify the regions where grid resolution is high enough to allow for hybrid 
simulation such as PANS. Also, we extract information about which regions of the flow we must set to f k 
= 1.0 in order to maintain a RANS simulation. We have come to the decision that the most appropriate 
parameter for assessing the grid resolution is the characteristic turbulence length-scale ratio. 

Here, we will highlight the basic concepts of this approach. Based on a simple dimension analysis, we 
assume that the turbulent viscosity may be related to the total kinetic energy (k),e, S and A as: 

v t « ^«A 2 S«A 2 — (8) 

£ K 

k 3 ' 2 
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L t is the characteristic turbulence length scale, which leads to 
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A 
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C h is a model coefficient of order one, which needs to be calibrated. In the present paper, we will use a 
value of one to evaluate the model. We will introduce two forms of the two-stage approach. First is the 
sequential form, which can be applied to wide range of unsteady flow conditions. We have successfully 
implemented and tested this form in the PAB3D CFD code. The second form is more complicated and 
requires major changes in CFD codes and needs larger resources. However, this second form offers a 
greater potential for solving more complicated flow problems such as moving bodies. 

A. Sequential Two-Stage Approach (PANS-v2) 

The simplest form of the two-stage procedure is shown in Figure 2a. The users will need to solve 
either the time-accurate or steady state RANS equations depending on the flow conditions. In the case of 
unsteady flow simulation, it is necessary to compute the time-average of the characteristic length-scale 
ratio. Equation (10) is then used to calculate f k . During the second stage (the hybrid method phase), this 
value will be fixed in time but it will vary in space. The users can quickly find out if the grid resolution is 
appropriate to carry out the PANS or other hybrid simulations. The second stage of this procedure is not 
in any way limited to use with PANS. In fact, any of the hybrid methods could use the results of the first 
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stage to construct a primary filter for use with these methods. This filter can be used to identify and 
isolate the RANS regions during the second stage of the simulation. In the present paper, we use the two- 
equation model described in the previous section to provide the first-stage solution. We refer to this 
procedure as PANS-v2. 

Here, we summarize the guidelines to be followed for the sequential two-stage procedure. These 
guidelines are completely dependent upon flow complexity. For the first stage: 

1 . Complete a three-dimensional, or two-dimensional simulation 

2. Unsteady or steady calculation, high-order schemes are not required 

3. Desired level of an allowable RANS turbulence model (one-equation, two-equation, Algebraic 
Stress, full Reynolds Stress, ...etc.) 

For the second stage: 

1 . Three-dimensional simulation 

2. Unsteady calculation, high-order schemes should be considered 

3. Hybrid models used (DES, Hybrid RANS/LES, PANS, . . .etc.) 

The users need to use the same flow conditions, boundary conditions, and grid resolution for both stages 
of the procedure. 

B. Concurrent Two-Stage Approach (PANS-v3) 

The concurrent form of the two-stage procedure is shown in Figure 2b. Users will need to solve RANS 
with standard turbulence model (STM) equations and time-averaged equations such as PANS with hybrid 
turbulence model (HTM) equations simultaneously during code execution. From the first-stage, we 
extract the instantaneous characteristic length-scale ratio and produce f k . During the second-stage (the 
hybrid method phase), this value will be used with the selected HTM to close the PANS equations. This 
process is repeated for each time step of the flow simulation. This approach is not demonstrated herein. 

Results and Discussions 

Flow past a stationary cylinder was computed at a Re =50,000 to validate the implementation of an 
efficient two-stage procedure for solving hybrid turbulence models. The grid consisted of 1,290,240 cells 
and 24 blocks, extended 15 diameters into the far field in the radial directions, and covered a two-cylinder 
diameter in the span-wise direction. This grid was developed and used by Vatsa et al. [14]. The same grid 
was used for all runs, which gave a first grid height y+ range of approximately .2 to 2.0 around the 
cylindrical surface. A non-dimensional time step of 0.015 (based on free stream speed and the diameter of 
the cylinder) was used. Based on the time step, approximately 350 frames of solutions in time per cycle of 
shedding were sampled. Four sub-iterations per time-step were used to reduce the error. The cases each 
required approximately 48 hours of CPU, using 24 (2.8 GHz P4) computers. 

The Laminar Separated (LS) case was chosen to evaluate the implementation of the PANS 
formulations into PAB3D code, and to investigate their capability to simulate such flow. References 
12,13 and 14 each gave a complete description of this case and the different computational procedures. 
This case is compared with the experimental data of Cantwell and Coles [15]. We followed the procedure 
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as described in Figure 2a. First, we used the RANS formulation to get time-averaged quantities to 
calculate the characteristic length-scale ratio (k as shown in figure 3a). This ratio varies in space and is 
used to produce the unresolved kinetic energy parameter (f k ). Figure 3b shows the distribution of this 
function at a selected spanwise (y) location. This parameter identifies the RANS and PANS regions. The 
RANS regions are defined with the parameter set at a value of one. The PANS regions are, for all the 
values, less than one. We use this parameter in solving the PANS formulation. Figure 4 shows the 
unresolved turbulent kinetic energy produced with the use of RANS (4a) and PANS (4b) formulations. 
The results of the RANS simulation were ten times those produced by the PANS formulation. More 
energy was resolved using the PANS formulation as compared with that resolved by the URANS 
simulations. 

Figure 5 shows the power spectral density (PSD) as a function of the Strouhal number, St, from RANS 
and PANS results. It is clear that RANS produced a peak PSD level at only one major Strouhal number 
and no significant minor ones. This result indicates that the 3D RANS simulation has a two-dimensional 
flow character. The vorticity plot shown in Figure 6 supports this conclusion as it shows the carpet-like 
flow with no change in the Y-direction. On the other hand, the PANS result shows PSD peaks at one 
major Strouhal number and several minor peaks at other St values, indicating that more scales were 
resolved in this simulation. Figure 7 shows the three-dimensional character of the flow produced as a 
result of the PANS formulation. 

The time-averaged surface pressure distributions from the 3-D simulations are compared with the 
experimental data of Cantwell and Coles [15] in Figure 8. Table 1 shows the comparison between the 
recent PAB3D simulations, the DES solutions of Travin et al. [13], and the experimental data. The present 
PANS formulations produced results comparable to the DES Travin et al. [13], and experimental data. 


Table 1. Re=50,000 LS Time- Averaged Results 


Method 

C D 

-Cpb 

St 

3D-DES Travin et al. [13] 

1.27 

1.28 

0.21 

2D-DES Travin et al. [13] 

1.77 

2.05 

0.14/0.2 

3D PAB3D PANS-vl 

1.1 

1.03 

0.21 

3D PAB3D PABS-v2 

1.14 

1.18 

0.21 

2D/3D PAB3D RANS 

1.08 

1.03 

0.23 

Experiment 

1.25 

1.2 

0.18-0.21 


Summary and Future Work 

We have introduced and implemented a novel two-stage procedure to efficiently estimate the 
level of scale resolution possible for a given flow on a given grid for Partial Averaged Navier-Stokes 
(PANS) and other hybrid models. The new capabilities improve the accuracy and robustness of creating a 
simulation of an unsteady flow field. In this paper, a new approach to prescribe the f k function was 
presented. This is a two-stage procedure, which is designed to remove the grid dependency of the hybrid 
RANS models. In the first- stage, we use the RANS simulation with a Standard Turbulence Model (STM) 
such as ke to produce an estimate of f k over the entire grid domain. In the second-stage, we supply f k for 
the selective application of a Hybrid Turbulence Model (HTM) such as PANS formulation in regions 
where the grid density is sufficient to resolve a portion or all of the large-scale flow structures. In the 
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present implementation, f k is a function of length-scale and grid size, which represent a characteristic 
length-scale ratio. For the simple two-stage approach, this function only varies with space throughout an 
unsteady flow simulation. 

In this technical memorandum, we selected the Laminar Separated (LS) case to examine the capability 
of the PANS formulation. With the use of the simple two-stage procedure an improved prediction 
capability is achieved in predicting C D , St, and the surface pressure distribution when compared to 
experimental data. This implementation is a first step towards adding a variable resolution turbulence 
model capability to CFD codes. We can now use the PAB3D code to refine the PANS formulation and 
to conduct validation computations using a variety of simple and complex flow physics problems. 

The present procedure represents the first step in the development and implementation of a more 
comprehensive approach for solving unsteady flows. This approach needs to be calibrated, verified and 
validated for a wide range of flow problems such as jet, cavity and others. Flows with moving bodies will 
present a great challenge for the present approach. This challenge could be addressed by extending the 
present procedure to the concurrent two-stage form. The concurrent two-stage procedure will require 
solving RANS with Standard Turbulence Model and time-averaged equations such as PANS with Hybrid 
Turbulence Model simultaneously at each time step for a total of 14 equations. This approach will need 
large resources of computers and code development for successful implementation. 
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Computation Domain = RANS + PANS Regions 



PANS Regions 


RANS Regions 


Figure 1. Regions for the use of PANS and RANS 




Figure 2a. Sequential Two-Stage Procedure for Solving Hybrid Turbulence Models 



Figure 2b. Concurrent Two-Stage Procedure for Solving Hybrid Turbulence Models 
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to support Hybrid 
RANS Simulation 

(PANS) 


Figure 3a. Characteristic Length-Scale Ratio (X), Contours at Y=1.0 
produced from RANS Simulation (Re=50,000) 


f k = 1 RANS Region 



f k < 1 PANS Region 

Figure 3b. Unresolved Turbulent Kinetic Energy Parameter, f k Contours at Y=1.0 Produced 
from RANS Simulation (Re=50,000) 
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Figure 4a. Unresolved Turbulent Kinetic Energy, k u (m 2 /sec 2 ), Contours at Y=1.0 Produced 
from RANS Simulation (Re=50,000) 



Figure 4b. Unresolved Turbulent Kinetic Energy, k u (m 2 /sec 2 ), Contours at Y=1.0 Produced 
from PANS-vl/v2 Simulation (Re=50,000) 
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St 

Figure 5. Power Spectral Density vs Strouhal Number Results from ke 
RANS and PANS Formulations (Re=50,000). 



Figure 6. Three-Dimensional Vorticity Magnitude Results from 
RANS Simulation (Re=50,000) 
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Figure 7. Three-Dimensional Vorticity Magnitude Results from 


PANS Formulation (Re=50,000). 



Figure 8. Coefficient of Pressure on the Cylinder Surface using ke RANS and PANS 
formulations Compared with Experimental Data 
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14. ABSTRACT 

The main objective of this article is to introduce and to show the implementation of a novel two-stage procedure to efficiently estimate the 
level of scale resolution possible for a given flow on a given grid for Partial Averaged Navier-Stokes (PANS) and other hybrid models. It has 
been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The first step is to solve the 
unsteady or steady Reynolds Averaged Navier-Stokes (URANS/RANS) equations. From this preprocessing step, the turbulence length-scale 
field is obtained. This is then used to compute the characteristic length-scale ratio between the turbulence scale and the grid spacing. Based 
on this ratio, we can assess the finest scale resolution that a given grid for a given flow can support. Along with other additional criteria, 
we are able to analytically identify the appropriate hybrid solver resolution for different regions of the flow. This procedure removes the grid 
dependency issue that affects the results produced by different hybrid procedures in solving unsteady flows. The formulation, 
implementation methodology, and validation example are presented. We implemented this capability in a production Computational Fluid 
Dynamics (CFD) code, PAB3D, for the simulation of unsteady flows. 
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